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Let M be a (i dimensional Riemannian manifold, with metric tensor ds"^ = 9iji^) 

l<i,j<d 

Q ■ dxi®dxj and consider a particle of unit mass moving along geodesies of M with kinetic energy 
^ I T = ■ Suppose we are now given a potential V{x). There is a relation between the 

geodesic flow and the Hamiltonian flow parametrized by t G R, for Hamiltonian H{p, x) = 
^ . 2 9^''{x)PiPj + V{x) on the energy surface E. To find the dependence {P{t), X{t)) on 

l<i,j<d 

the new time t one has to solve the equation dt = {2{E — V))~^^'^ds which we call Maupertuis 
principle (see [Ar2]). Then the particle moves along the geodesic flow on M relative to the 
metric dp = {E—Vy^'^ds, and we say that both flows (and metrics) are conformally equivalent. 

We shall consider in this paper the following stronger variant. Maupertuis Jacobi cor- 
respondance for the pair {H.Ti), T-L{p,x) = 2{e-v) 9^''{^)PiPj: ensuring that both 

l<i,j<d 

systems have the same trajectories {P{t), X (t)) and {X{t),V{t)) at energies {E, 1), is given 
by dr = [E — V)dt. This principle was used for instance by T.Levi-Cevita in relation with 
Kepler problem, cf. [Ts], [HieGrDoRa] and references therein ; see also [EasMat] for the 
case of projective equivalence of metric connexions. Among interesting situations of applica- 
tion of Maupertuis-Jacobi correspondance is the case, when the functions (A'(r),P(r)) are 
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T-quasiperiodic : {V{t), X{t)) = {J^{ujt + <^°)), X{ijjt + <^°)). Here a; is a vector of periods 
and (7'(t), A'(t)) are 27r-periodic functions with respect to each phase. 

This occurs when V. is completely integrable, and A are compact, connected leaves of 
maximal dimension d (i.e. lagrangian) in the fibration of T*M by the momentum map m, 
that we assume to be non singular at A = Aq. Then Arnold- Mineur-Liouville theorem states 
[Ar2] that in a suitable action-angle coordinate system, A are diffeomorphic to the torus 
rpd _ R,'^/27rZ'', and parametrized by the action variables / e R^. The non-criticality of m 
implies that all components Ij of / are non zero in a neighborhood of Aq. 

In general the new hamiltonian H has no reason to be integrable. Nevertheless, in a 
neighborhood of Aq, there are action-angle coordinates (J, (p) giving a full family of tori A"^ 
which are almost invariant for the hamiltonian vector field Xh- Under some diophantine con- 
dition, it it possible to improve these tori, so that they become invariant within an arbitrary 
accuracy. 

This is actually the case, for using a generalization of Birkhoff transformations (better 
known near a stationary point of Xh, ) in the symplectic coordinates given by Darboux- 
Weinstein theorem, H can be brought to a normal form, we call the Birkhoff normal form 
(BNF). This was carried out by M.Hitrik, J.Sjostrand and S.Vu-Ngoc [HiSjVu]. 

An alternative, and more pragmatic way to our opinion, was devised in [BeDoMa], and 
consists in providing directly suitable action-angle coordinates (J, without resorting to 
Darboux-Weinstein theorem. We believe it has the advantage of being easily implemented 
numerically in some concrete applications. 

Once the classical Hamiltonian is taken to BNF, we may address the problem of finding 
quasi-modes supported on the A"', with a suitable accuracy. A necessary condition for A'^ is 
to satisfy Maslov quantization condition. A similar problem arises in an attempt of quantiz- 
ing the so-called KAM tori. Quoting J.J. Duistermaat [Dul,p.231], while slightly changing 
notations : 

"... there is no reason why the tori on which [Maslov quantization condition] holds for 
some h > should persist. For most [invariant embedded lagrangian manifolds i : A — T*M] 
one can expect that h ^ j^- 6 [the cohomology class of the closed 1-form i*{pdx)] is dense in 
Jf^(A; R)/i7^(A; Z). Extracting a subsequence |- = r^^,, from = tq +ncro,n G Z [which 
ensures Maslov quantization condition to hold], such that ■ 9 converges to —a/4 [Maslov 
class] in H^{K; R)/if^(A; Z), one can still obtain some asymptotic estimates for the spectrum 
of H but in general the convergence of r^-^ to — q;/4 is expected to be slow and the asymptotic 
estimates for the spectrum are correspondingly weak." 

This observation was made more precise by Y.Colin de Verdiere [CdV,Prop.7.5] : namely, 
if i?^(A; R) is of dimension 2 (i.e. for a 2-d torus A), then the extraction above can be made. 
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provided ^, expressed in a basis of entire classes, is irrational. 

A different point of vue consists in selecting among the family A-^ suitable tori on which 
Maslov quantization condition holds, i.e. a sequence J = Jk{h) , k E Z, \k\h < . This allows 
to produce an almost "full" (or "density 1" ) spectral sequence for H, in some suitable sense. 

The first route, followed in [HiSjVu], consists, through /i-Pseudo Differential Calculus, 
in quantizing the BNF microlocally near the A"^, and this can be achieved for corrections to 
all orders in h. (See also [CdVl] , [Laz] , [Po] , and [Vu] for the case of a family of d commuting 
/i-PDO. ) 

This is a rather remarkable result, since quantization of manifolds is in general impossible 
beyond first order. It relies essentially on the fact, observed a long time ago by Colin de 
Verdiere [CdV], and used later by G.Popov [Po], that the action of H on semi-densities on tori 
verifying Maslov quantization condition, reduces to the action of another /i-PDO, obtained 
by conjugating H hy a Fourier Integral Operator (FIO) quantizing the BNF, to microlocally 
defined functions on the flat torus T''. These functions satisfy Floquet periodicity condition : 

(0.3) u{x -u) = e2^^<'^+"'^/^'">/'^«(a;), u e 

Here / = (/i, ■ ■ ■ , Id) are the classical actions parametrizing the torus A (stated otherwise, / 
is a rcprcscntant of Liouville class 9 in H^{A; R)), so that u{x — u) — u{x) for all integer u 
whenever 9/h + a/4 = mod Z^. 

Here we have chosen another route, consisting simply in using the almost invariant 
tori A"' to construct quasi-modes for operator H{x, hD^) near energy 0. This can be done 
by adjusting the new action variable J so to satisfy Maslov quantization condition, and 
solving the transport equations for H to first order along A"'. We believe again that these 
constructions are simpler from a pragmatic point of view, especially when H assumes the form 
H = —h?/S.-^V. In more complex situations however where (see Section 4.d) we have to resort 
to FIO's implementing the change of variables that takes our operator to the normal form, 
passing for instance from a general conformal metric to a model one, but the constructions 
are very similar. 

We could also include isotropic tori, corresponding to the gaps in the foliation of the 
energy surface by lagrangean manifolds when the system is no longer integrable, or consider 
critical points of the moment map, which leads to introduce so-called harmonic oscillator 
coordinates (in the elliptic case). 

The paper goes as follows. In Sect.l we recall from [BeDoMa] the construction of action- 
angle coordinates in the neighborhood of a torus A invariant for the hamiltonian flow Xh, 
leading to existence of a family of tori A'^ '"^ in the neighborhood of A which, if A is diophan- 
tine, are quasi-invariant over arbitrarily long times. Thus we get a variant of BNF. 
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In Sect. 2 we quantize some of these A"^ '^, 



so that they satisfy Bohr-Sommerfeld-Maslov 



quantization condition, and support quasi-modes for H with a suitable accuracy. 

In Sect. 3 we apply these considerations to semi-classical Jacobi-Maupertuis correspon- 
dance : assuming that H. is completely integrable, or simply has an invariant torus A in some 
energy surface, we show how to construct quasi-modes for H. We show that the frequencies 
uj (resp. a;) of the quasi-invariant flows for 7i (resp. H) are simply related by a; = {g)oj. This 
implies partial isospectrality of both hamiltonians. 

In Sect. 4 we give various applications of construction of action-angle variables, and exam- 
ples for semi-classical Jacobi-Maupertuis correspondance, including the shallow-water waves 
theory. Namely we show how to select profiles of the basin, corresponding to a family of 
Liouville 2-D integrable metrics, which give raise to trapped modes of a given wavelength. 
This could be easily implemented on a computer program, and consists actually in the main 
"practical" result of that paper. 

In Appendix, we give a short account for Maslov theory and WKB constructions as 
needed in the main text. 

Acknowledgments: One of us (M.R.) thanks Ch.Duval, J.Sjostrand and D.Hafner for useful 
remarks. We also acknowledge a grant from the Region Provence- Alpes-Cote d'Azur in 2007, 
which enabled to spend some fruitful and enjoyable time in both institutions. 

1. Action-angle variables for a family of almost invariant tori. 

First we recall the main result of [BeDoMa] . Consider a smooth Lagrangian torus Aq = 



{{p,x) = (P°((^),XO((^)), ip= ((pi,---,(pd) e T*^}, = (R/27rZ)'' embedded in T*R'', and 



be the corresponding actions computed along the basis of cycles (7^ , ■ ■ ■ , 7^) on Aq. We denote 
by Y^{<f) the 2d vector function *(P'^, X^). Although it is possible to take degeneracies into 
account, we shall assume here that the moment map is non singular near Aq, i.e. /? 7^ for 



We may also assume that Aq belongs to a ^-parameter family of tori A*, s varying in an 
open set of R^, 1 < i < d. This is the case, by Arnold-Liouville-Mincur theorem, if these 
tori arise from an integrable, or semi-integrable Hamiltonian system. If £ = 1 we can think 
of s = E as the energy of the system, but more generally, the family can be parametrized 
(possibly after renumbering variables, ) by the first £ action variables {Ii, ■ ■ ■ , Ii). However, 



let/° = (/°,---,/°). 



(1.1) 




all J. 
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tori of the family A* are not necessarily preserved by the transformations below, so to avoid 
technical diflficulties we shall here take simply £ = 0, and consider a single torus Aq. 

Introduce the 2d x d matrix with columns {Y^^ , • • • , Y^J, and the d x d symmetric 
matrix M = —{^Y^Y^)~^ (we denote often partial derivatives by a subscript. ) Let also J 

( -Idd\ 

be the standard symplectic matrix J — \ , where Id^ is the d x d unit matrix. 

Vidd ) 

We put 

(1.2) = = JY«M 

The 2d X d matrices Y^ and Y"^ satisfy the conditions 

(1.4) 'Y^JY^ = 0, ^Y-^ JY^ = -Idd, 'Y'^JY-^ = 0, 

the last two being implied by (1.2) and the first one, which precisely express that Aq is a 
Lagrangian torus. This is proved in Appendix of [BeDoMa] by a direct computation. 

In the notation Y"^ , index J refers to new action variables we are constructing below, for 
which Y-^ enjoys the same symplectic properties as ^ for an integrable system. All these 
matrices are periodic with respect to (/? e T''. 

Now we introduce precisely the new action-angle coordinates (J, </?) G neigh(/'^, R'^) x 
(R/27rZ)'^, so that y' = Y{J,(p) — (x(j^)) canonical coordinates defined in a neighbor- 
hood of Aq. As is shown in [BeDoMa] this can be achieved by taking 

(1.5) Y{J,ip) = Y'{ip) + Y'{^)G{J,^) 

where the real d- vector G — *(Gi, • • • , Gd) is obtained by solving the system of the quadratic 
equations : 

(1.6) ,- Gj - ^{G, 'Y'JY^.G) = 6„ J = 1, 2, . . . , d 

with data l = (ti, ■ ■ ■ , td) = J — I^. We fix the root of this equation by the requirement 
G{J,(f) = L + C(t^) as i — > 0. Equations (1.6) can be solved with the help of a fixed 
point theorem, and G depends smoothly and periodically in (p, analytically in J (or t). For 
completeness we recall the following : 

Lemma 1.1 [BeDoMa]: Coordinates y' = {p' , x') = y + 0{l) are symplectic, i.e. they verify 

(1.7) (y;,,Jy;j = (yj,,JyjJ = 0, {y'j.,Jy'^,)^-Sjk, j = i,...,d 
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Proof. Consider first {y'j.^Jy'j^)- Differentiating (1.5) w.r.t. J and making use of the last 
relation (1.4) we find {y'j.^y'l) = {Gj-.^Y^JY^Gj,) = 0. 

Consider next Ajj. = {y'j., Jy'ip^) ■ Differentiating (1.5) w.r.t. (p^ and Jj, then making 
use of the second relation (1.4) in the form ^Y^ JY^^ = — Cfe, being the A;:th unit vector of 
R'', and of the last relation (1.4), we find 

Ajk = -Gkj, + {Gj.,'Y'jY^^G) 

Then differentiating (1.6)j w.r.t. J^, and making use of the relation ^Y^ JY^, — —^Y^.JY"^ 
which follows from the last equation (1.4), gives Aj^ = ~^jk as stated. 
At last, consider B^j — {y'^i^, Jjy'^_.), which by (1.5) again gives 

(1.8) Bkj = {JY^. +JY^.G + JY-'G^^,y;^+Y^^G + Y-'G^,) 

Differentiating (1.6)^ w.r.t. (fkj then (1.6)^ w.r.t. to (fj, and forming the difference, the 
mixed derivatives Y^^^. cancel out, which gives : 

~^jVfc + ^k^fij + 2^G^^,*Y-^JY^.G) - -{G^.,*Y-^JY^^G)+ 

1{G, i^Yi^JYi, - 'Yi,JYi^)G) + \{G, 'Y^J^Yip^, - F/^G^J) = 

We have 'Y^^JY^. - 'Y^.JY^^ = S^^i^^i- {G^'Y^JY^.G^^) = (G^, , ^F^JTF/.G), 
so 

-G.v, +G',^, + {G^jY'^JY^p) - {G^^,'Y'jY^p) + (G, *y/^ JF/^.G) = 

Replacing the term (G, ^Y^^JY^.G) in (1.8) by its expression, we see that {G^. , ^Y-^ JY^JG) 
cancels out with {^Y^^JY^ G^.,G), {G^^.^Y^ JY^.G) cancels out with {^Y^ JY^.G,G^^) , 
while {'Y^JY^G^.^G^^) = {JY^.,Y^J = by (1.4). So (1.8) simplifies into 

Bkj = Gj^p^ - Gk^. + {JY^.G, F^J + {JY-^G^. , + {^Y^^JY^. , G) + {^Y-^ JY^. , G<pJ 

and because again of ^Y"^ JY^^ — — e/j, B^j = 0. The Lemma is proved. J(t. 

In other words, the transformation k : {J,f) ^ y' = K{J,ip) is canonical. This yields a 
(i-parameter family of tori A"', by choosing J = + l— Const., and A"' = Aq + 0{l) in the 
C°° topology. 

Remark: Actually wc know, by Darboux-Wcinstein theorem (sec [BaWe, Corollary 4.20]) that 
if A C T*R'^ is a lagrangian submanifold, then a tubular neighborhood of A in T*R'^ is 
symplectomorphic to a neighborhood of the zero section Z a; A in T*A, by a map which is 
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the identity on A. The construction above speciahzes and provides a simple proof to Darboux- 
Weinstein theorem in the case where A is a torus given by J = /° in action-angle coordinates, 
so that neighborhoods of Z in T*A can be identified with neigh(/°; R^) x T'^. In other words, 
given the lagrangian immersion i : A — > T*M, we have constructed a symplectic immersion 
X of a neighborhood of Z into T*M, such that T = i o tta on Z, where tta : T*A ^ A is the 
natural projection. 

We still denote by the basic cycles on A'^, which are obtained by a continuous defor- 
mation of those on A, because of the nondegeneracy condition. So we have 

(1.10) Jj = — l p'dx' 

Assume now that Aq is an integral manifold for Hamiltonian H{p, x), i.e. p = P°(<^ + a;t), x = 
X^{(p + ujt) is a quasi-periodic solution of the hamiltonian system (p, i) = Xh{Pix) on Aq 
with rationally independent frequencies a;, we try to construct a normal form for H near 
Aq. Here notation Xh{p,x) stands for the Hamiltonian vector field of H{p,x). By a Taylor 
expansion in coordinates y' as above, we have : 
(1.11) 

H o k{J, ^) = H{Y\^) + y^(<^)G'( J, ^)) = H{Y\^)) + {H'{Y%^)), {^)G{J, ^)) + 0{i'') 
with H'(Y^{ip)) = -Jy = -J{u;, i^)Y\^). So 

{H'{Y\^)\Y\^)G{JM) = Y.u>,Yf{^)Yi{^)Gm{JM 

which by (1.4) and the relation G(J, ip) = l + 0{i^) gives 

(1.12) H{y') = Ho + {uj, i) + 0(^2), Ho = H\^^ 

so the A"' are invariant at first order in i under the fiow oi H o k. 

Assuming further a diophantine condition on w, we can extend this normal form as 
follows, using a finite order KAM procedure. We look for a generating function S{J',(f) 
such that J = d^S{J',ip), ip' = dj>S{J',ip), where J' = + i! . We find that S{J\^) = 
(J',(/?) + $(/-',(/?) is defined as a Fourier series in G T'^, with vanishing zero mode, and 
coefficients $fc(t') = C(i'2). Push Taylor expansion (1.11) to next order in t = t' + ^(^$(6', (/?), 
we find after averaging over Lp the quadratic terms in l' 

(1.13) H{P{J, cp),X{J, cp)) = Ho + (a;, l) + A{l') + B{l', cp) + {u, d^^t', cp)) + 0{l'^) 

where A{l'), B{l' , (p) = 0{l''^) are given smooth functions defined in a neighborhood of Aq, 
and {B{l', •)) = (average over T*^. ) Since u> is diophantine, and B{l', (p) smooth in (p, we 
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can determine the quadratic part of $ in i! by solving B{i'^(p) + {uj,d^^{L' ,(p)) = up to 
0{l'^) terms. This process easily carries out at any given order, yielding the following : 

Theorem 1.2 [BcDoMa]: Let H be as above, and assume the frequency vector uj satisfies 
a diophantinc condition. Then for any iV = 1,2, ■ ■ -, there is a smooth canonical map kn : 
(J', (/?') (J, Lp) defined for J' in a neighborhood of / and ^p' G T'^, and a polynomial Hn{l') 
of degree N, with HN{i') ^ Hq + {u, i') + 0{i''^), Hq = (w, /) such that 

(1.15) HoKo kn{J', (fi') = HN{i') + 0{i'^+'^) 

This is precisely a Birkhoff normal form for H in o. neighborhood of Aq (see also 
[HiSjVu]. ) Assuming analyticity on these estimates could certainly be sharpened to 
give a remainder 0{e~''' ^/^) as in Nekhoroshev theorem, < /? < 1 being related to the 
diophantine condition, and C > (see e.g. [Po2]. ) 

To conclude this Section, we make a few remarks concerning quantization in action-angle 
variables, following [BeDoMa]. Recall [BaWe, Definition 5.33] that a linear map p from the 
space C°°(T*M) of smooth functions, to the algebra A generated by self-adjoint operators 
on some complex Hilbert space and endowed with the Lie algebra structure defined by 
[A, B\h= ji{AB — BA), is called a quantization provided it satisfies so-called Dirac axioms : 
(1) p(l) = Id, (2) pi{f,g}) = [p(f),p(g)]h, (3) for some complete set of functions /i,- • •,/„ 
in involution, the operators p(/i), • • • , p(/n) form a complete commuting set. We know this 
set of axioms is in general too stringent, in the sense that a quantization of all classical 
observables doesn't exist, although the symbolic calculus of /i-Fourier integral operators, or 
Maslov canonical operators, provides a fairly good approximation of this classical-quantum 
correspondence. 

Quantization deformation occurs already in the simple case of a completely integrable 
hamiltonian system on T*M, with hamiltonian H{p, x), which admits a family of lagrangian 
tori A^. Namely, trying to quantize the corresponding action-angle variables (/, (p), considered 
as classical observables, in a neighborhood of the A^'s, we require that p{I) = I.,p{(p) = (p 
would satisfy 

(1-17) [/j,/fc] = 0, [(pj,Ik]h = Sjk-, [(pj,(pk]h = ^ 

and moreover, that the semi-classical hamiltonian H{x,hDx) associated with H[x,p) via 
usual Weyl /i-quantization, would be a function of {(p, I). The naive answer would consist in 
choosing (p as multiplication by and / = hD^p. But then if we try to recover the canonical 
operators Xj = Xj{(p, I),pk = Pk{^,I) by symbolic calculus, it turns out that the canonical 
commutation relations \pk,Xj]h = Sjk are only satisfied modulo 0{h). (This can be easily 
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checked also when H is the harmonic oscillator, and x = Vlcos(p,x = y/l sin (p, H{p, x) = 
+ x^ = I. ) Therefore one should in general also include quantum corrections and seek 
for canonical variables as Xj = Xj{(fi,I) + hXj^\(fi, I) + • • •. This in turn affects the full 
symbol of hamiltonian H, expressed in action-angle variables, by a quantity 0{h^), and 
H{p,x,h) = Ho{I) + h'^Hi{I,(p) + ■■• depends again on angle- variables. But introducing 
e = as another small parameter, we can apply a KAM procedure, under some non- 
resonance condition, to replace (/, (f) by new symplectic coordinates (/', ip') and thus get rid 

of angle dependance in the 0{h^) term, which yields H{p, x, h) = Hq{I) + h?'Hi{I) H . See 

[CdV2], [Vu] for the case of commuting /i-PDO's. 

For a quasi-integrable system, we can combine the BNF (at the level of principal symbol) 
with the procedure above. 

In our case however, mostly because we content with the principal and subprincipal 
symbol of we have no real need to recover the canonical variables (p, x) from (/,</?), or 
(J, (/?), nor to quantize action-angle variables. 

2. WKB solutions along perturbed tori. 

In this section, we carry out WKB constructions associated to the family of perturbed 
tori A"^. We start from tori A-^ on which H verifies (1.12), and then extend the argument for 
tori obtained as in Theorem 1.2, which are closer to be invariant by the hamiltonian flow. 
We use essentially known results on geometry of tori that we have recalled in Appendix. 

a) Action integrals. 

We first compute the action integrals on the lagrangian tori A'' C T*M, constructed 
from A with J = / + z,, in terms of the action-angle variables (J, (p). We denote by JF the set 
of focal points, that projects onto the caustics C (see Appendix A.) 

To start with, consider a, 6 G M be away from C. Assume y = {p{b),b) G A is not a 
focal point, i.e. ^^|(^=(p/(6) is not singular ; here f^{b) is the angle-coordinate of b on A. For 
J G neig h(/), let Lj : T'^ ^ T*M be the la grangian embedding with image A"' (the context 
prevents from confusing the map Lj with the action variables z,), dc/)"^ = L*j(pdx) the 1-form 
parametrizing A"^ near 6, and Lp"^ (x) the angle-coordinate of x G neigh(6) on A'^. 

Let 7 C A have its x-projection contained in the domain of definition of the (j)"^ . We want 
to compare f^pdx with f^,p'dx' on 7' C A'^, where y' = {p',x') = Y{ J,(p) are symplectic 
coordinates as in (1.5), the x-projection of these paths joining a to b. So 7 and 7' are 
parametrized by 99, the action variable being held constant. Along 7, x = X^(Lp),p = P^(lp), 
while along 7', x' = X^{ip) + X-^ {(p)L + 0{i^), and p' = P^{ip) + P-^ {(p)i + 0{i'^), with A:'^(^) 
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and P-^{(p) as in (1.2). Substituting these values into J ,p'dx' we get 

/ p'dx'= f {P\^),dX%^))+ f d{P\^\X\^)i) 

+ / {P\^)i,dX\^))-{dP\^),X\^)i) + 0{i^) 

We choose the origin of angular coordinates such that (p'^{a) = 0. Because A"^ is lagrangian, 
Jyp'dx' only depends on the end-points of 7'. We have 



(2.2) / d(P°(^),X^(^).) = (P°(^),X-^(^).)|o^ 

Next we want to replace the path 7' in J^,{P^{(f),dX^[(p)) by 7, which we complete into a 
1-cycle in the (/9-space as 7' = 7'" + 7 + 7", with 7'" between and <f^{a), 7" between (f^{b) 
and (f^ib). By Stokes' formula and since dP^{(p) A dX^{(p) = : 

(2.3) [ {P\ip),dX%ip))^{[ + [ +[ ){P\^),dX\^)) 



We express the fact that b = X(^-^(6)) = X^{ip\b)) as 

(2.4) ( -^{^{s))ds,^\b)-^\b)) = -X\^\b))G{J,^\b)) 

where (/7(s) = (p^ {b)-\-s{(fi'^ {b)~(fi\b)) is a straight line. By continuity, <~p{s) 'p\b) uniformly 
on s G [0, 1] as t — 0. Since b is not a focal point on A, ^^{(f^ (b)) is not singular, and the 
same holds for ^^{ip{s))ds when l is small enough, so (2.4) gives ^'^{b) — ^^{b) = 0{i). 
Using also ^{^{s))ds - ^{ifi{t)) = C(t) uniformly on t G [0, 1] and P^{ip{t))dt - 
pO(<^^(6)) = 0{i) we find 

/ {P\v),dX\ip)) + (P0(^),X'^(^)6)|^=^.(,) = 0(6^) 

The same holds for the contribution of 7'" since we may always assume, a being fixed, that 
it is not a focal point, so these estimates shows with (2.2) that 

(2.5) f d{P\^\X\^)i) + {l +[ ){P\^),dX\^)) = 0{t') 



On the other hand, using the canonical relations (1.7) for coordinates y', and the estimate 
G{J, (f) = i + O(i^) in C^-topology, one readily verifies that 

{P'{ip)i, dX\^)) - {dP\^), X'{cp)t) = {t + 0{l^), d^) 
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Integrating this 1-form along 7', where l is constant, and inserting in (2.1) together with (2.5) 
we have : 

(2.6) f p'dx' = f pdx+ {if\h) - <p\a),i) + 0{i'^) 

J ^' J ^ 

More generally, for N = 2, 3 • • •, kjv as in Theorem 1.2 and J' e neigh(/), consider the 
lagrangian embedding lj'^n = io kn : T'' — > T*M, with image A"^''^ . Let again dcf)'^''^ = 
i*j, j^{pdx) parametrize A'^'--'^ near h. Using the estimates (p' = (p -\- 0(i'^), J = J' + 0{i''^), 
one can readily extend (2.6) to the tori A"^ , so we get : 

Proposition 2.1: Let a be away from the caustics of A and y = {p{h),h) e A is not a focal 
point, and assume a and h can be joined by a path 7 C A contained in the domain where A 
is projectable (i.e. the x-projection of 7 is contained in the domain of definition of (^^, where 
the 1-form dcf)^ = i}{pdx) is a local parametrisation of A. ) Then for any N = 1, 2, • • •, and t' 
small enough, for all tori A"^''^ as above, we have : 

(2.7) 0-^''^(6) - 0-^''^(a) = 0^(6) - 0^(a) + (<^-^''^(6) - <^-^''^(a), + 0{l") 

Note that (2.7) can be written at the infinitesimal level as t}/ j^{pdx) — i}{pdx) + L'difi + 

Next we consider the case where y = {p{b),b) G A is a focal point, but y ^ J-'o [which 
we recall from Appendix A], so that ^^£^\ip=ipi{b)j (i-^- diTp : A — > R^) is non singular. 
For small z,, the "alternate" action is now f~ —xdp and the new generating function dip^ = 
i*j{—xdp). At points y where rfyr^ and dn^ are non singular (i.e. where A^ is bi-projectable) 
V'^ and are related by Legendre transformation (j)^ [x] = {—p^V'^^{p)) + i^\p)\p=p{x)^ so 
= •0^ + i*j[xp)^ up to a constant. We pick again a away from the caustics, and write 
y = (/9,a;(/9)), (p(a),a) = {a,x{a)). For sufficiently small t this holds for the neighborhing 
tori A'^ and also for the A'^ . 

So we need just repeat the former argument, comparing j~, —x'dp' with j~ —xdp, where 
the first integral is computed along a path 7' joining (the lift of) a with (3 on A"^, and the 
second one along a path 7 joining a with (3 on A, and contained in the domain of definition 
of ifj^ . Thus we get : 

Proposition 2.2: Assume y — {p{h),h) = {/3,x{/3)) G A is a focal point, but y ^ J^q. If 
a and b can be joined by a path contained in the domain of definition of i/j^ , then for any 
N = 1,2,- ■ ■, and t' small enough, for all tori A"' '-'^ as above we have : 

(2.8) - V'^{a) = ij\/3) - ij\a) + (^^''^(6) - V'^''^(a), t') + 0{i'^) 
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with notations as in Proposition 2.1, and where dip 
ametrizing A"^ near b. 



J',N _ 



tj, ^(—xdp) is the 1-form par- 



Note also that (2.8) can be written at the infinitesimal level as pf{—xdp) — L}{—xdp) + 



Then we consider the general case, to take in account the possible exceptional set JFq 
defined in Appendix A so we assume that dn\A{y) is of rank k < d — 1. So A can be 
parametrized by a Morse function S{x,6) defined near ij^{y), and by the discussion leading 
to (b.27), we may assume S{x, 6) is of the form T(xi,p2) + (2^25^2) and A = A^. 

Take a base point A = (02,01) G R'^ x R^~'^ that we identify with its lift through 
TT : A — R'^, (p, x) 1-^ (p2, 2:1) with A G At, we assume that A is non singular. The variable 
point near y, with the same abuse of notations is now B = (/32, bi) G At- We deform A into 
A"^, and the generating function T to T'^ . We shall always work in the domains where T'^ 
are well defined, and drop the subscript T"^ from A;^j. So let 7 C A be a path between (the 
lift of) A and -B, 7 = 71 +72 such that p2 = /?2 along 71, and xi = ai along 72. Using (b.27). 



we find f^^Pidxi = -§^{xi,P2)dxi\p^=(3^ and /^^ -X2dp2 = ^{xi,p2)dp2\xi=ai, so 



In the same way, we let 7' C A'^ be a path between A and S, 7' = 7^ + 72 such that p2 = P2 
along 7^, and xi — ai along 72, so denoting T' for T"^ we have : 



We parametrize A^ and A'^ by 99, so that x = {xi,X2) = X^{lp) = {Xi{(p), (X2((/7)), along 
7, and x' = {x'ijX2) = X{ip) = (Xi (</?), (X2(v')), along 7', where as above, X{lp) = X^{lp) + 
X"^ {Lp)i+0{i'^), and -P(v) = P^{(p)+P'^ {^)i+0{i'^). Possibly after renumbering the variables 
we may assume that we have a partition of the coordinates = ((/?i, (p2) adapted to the 
projection tt : A ^ R x R , i.e. such that and are non singular. Then the 
argument of Propositions 2.1 and 2.2 carries over to the present situation, working separately 
on 71 and 7^ with the xi and x'l coordinates, then on 72 and 72 with thep2 andp2 coordinates. 
So we proved the following generalization of Propositions 2.1 and 2.2 : 

Proposition 2.3: Assume y = {p{b),b) G A is a general focal point, let A = A-'^ is 
parametrized as in (b.27), and pick a non singular poiint ^ G A near y. Then for t small 
enough and B near y, 



(2.12) {T'{x[,p'2) + (4,P2))Ia = {T'{x^,p2) + {x2,P2))\^ + {v\b) - ^\a),i) + 0(6^) 



i'd^ + 0{i''^). 



(2.10) 




(2.11) 
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Here d{T-^ {x[,P2) + (3^2)^2)) = i^j{~xdp) = i^jipdx) is the 1-form parametrizing Aj near B 
where it is bi-projectable, and (p'^{B) e T'^ the angle coordinate of B. More generally, let 
N = 2,3, • ■ •, and l' small enough, for all tori A"^ as above we have with obvious notations : 
(2.13) 

(T-^''^(a;;,p^) + (4,P2))Ia = {T'ix^,p2) + (0:2,^2)) |f + ^(6) - <^^''^(a),0 + 
b) WKB solutions. 

The constructions are recalled for completeness in Appendix A in case of a Schrodinger 
operator —h'^A + V{x), but the result holds for more general /i-PDO. The local WKB solution 
near point € M depends again on the structure of A above xq. 

First we describe locally the quasi-mode associated with A^ in term of {P^ (if) , ((p)) . 
In the simplest case described in Proposition 2.1, i.e. where A is projectable, the phase 
function restricted to A equals (f)^ {(p), it verifies -P^(v') = ^4>^ {'P) ^-nd solves Hamilton-Jacobi 
equation H {V (l)\{p) , [ip)) = Hq = E. Recall from (b.4) that the solution to the transport 
equation verifies ao{x; -E) |x=xo(v3) = Const. | det ^^l""*^^^- 

Denote for short by pdx a solution of Hamilton-Jacobi equation with data on A, 
and let u{x;E,h) = ao{x; E)e'^r ^'^'"'^^^ . Then u solves {H{x,hD^) - E)u = 0{h'^). Here 
H{x, hDx) denotes any "reasonable" /i-quantization of H{p, x), e.g. Weyl quantization ; note 
that as usual, we first quantize the variables {p,x), then restrict to the torus. Moreover 
u{x; E, /i)U=xo(^) = Const. | det ^ I" V2e^'^^(x°M)/\ 

For A?" = 2, 3 • • • let us then consider, for small t' 

(2.15) u{x;c\h)U^^^j,,^,^=Const.\det^^^^^ ^' e 

Clearly, u{x; i' , /i)|x=x(j',v5') is a lagrangian distribution supported on A"^ , and 

\B.{x, hDa;)u{x; i' , h)] \x=x{j',^') 

is obtained by restricting H{x, hD^;) [ao{x; E)e'^I pdx)/h-^ tox = X{J' , cp'), andp = P{J' , (p'). 
Let also En = Hjqii')- At the level of the eikonal equation we have by (1.15) 

(2.16) {H{p,x) - EN)\y=^P^j',^'),xiJ',^')) = 0{i"'+') 

Now we consider the transport equation and compute, according to (b.3), (b.4) 

(2-17) Cx,\^^,Jao{x,E)\dx\'/')^Cx,^^,,^J\dv^\'/') 
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Here Cxh denotes the Lie derivative of Hamilton vector field Xh acting on half-densities, 
Xh is no longer tangent to A-^''^, but it follows again from (1.15) that Cx^ij^j j^{\d(f\^^^) = 
0{l'^+'^). So if we choose N such that t'^+'^ =0{h'^),we get 

(2.18) [{H{x, hD^) - En)u{x; l' , h)] U=x(j',cp') = 0{h'') 

Now we turn to the case described in Proposition 2.2, i.e. near a focal point yo = {po, xq) 
where A is TTp-projcctable. The phase function restricted to A equals i/;^ (^p), it verifies X^{ip) — 
V'0^(<^) and solves Hamilton-Jacobi equation i/(V'i/'^(v), -P°(v)) = Hq = E. The quasi- 
mode restricted to A near xq is then given by (b.5), with phase i/j^ , and principal symbol 
bo{p, E) = C\ det ^r^/^ that verifies (b.l5). 

So for N = 2, 3, • • •, and i' small enough, let x ^ is a cut-off equal to 1 near the 
angle-coordinate of po, substitute as above {P{J', ^fi'), X{J', (p')) for (p, x) in (b.5), and define 
(2.19) 

vix;c',h)\^=xiJ',,n = i^^h)-'^'J^/^^^'''^ 

Following the proof of Theorem b.l leading to (b.l5) we see that v{x; t' , h)\x=x{j' ,ip') also 
verifies (2.18), with the same energy En = H^ii') that depends only on A-^''^ . 

At last, we can consider the general focal point yo described in Proposition 2.3, and 
use the general representation of the quasi-mode supported on A in a neighborhood of xq 
given by (b.l7), i.e. w{x; E, h) = {27rh)~^/^ J e'^^^^'^^/^a{x, 6, E, h)dx, where we have chosen 
the number of phase variables equal to d. It is suitable to express w{x; E, h) in the mixed 
representation described after Theorem b.l, so take {xi,p2) as local coordinates on A, = 
ip{x\,p2)- For the unperturbed torus A, we have 

(2.20) X exp[i( |'(P,dX) - (P2(xi,p2),X,(a;i,p2)))A]|,=(P(,^,^^),^(,^,,^)) 

xx(^(x.P.))|det ^(^^y(^)V ^/^cfp. 

Then we argue as before, replacing (/, v?) by { J',ip') in (2.20), and verify it satisfies (2.18) 
when = Oih?). To get a quasi-mode we need to glue different lagrangian distributions 

expressed locally as tt, v, w and take care of Maslov indices as is explained in Appendix 
C. There we notice that the embedding (A, l) is topologically stable if we exclude some 
exceptional lagrangian singularities in A, and so Maslov indices for the A"^ remain the 
same. 
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On the other hand, Propositions 2.1-3 and imphcit functions theorem show that, given 
5 > 0, it is indeed possible to "match" i' with values of the actions and Maslov indices a along 
the fondamental cycles as in (1.10) so to satisfy Maslov quantization condition |- • J' + Q!/4 = 
mod iy^(A;Z) = H^{A.^ '■'^;Z), for a sequence J' = J']^{h). So let us summarize our main 
result in the 

Theorem 2.4: Let (A, i) as before such that A is invariant under the Hamiltonian flow 
of H and this flow has diophantine frequencies u, H\jy = Hq = E. Assume A is free 
of exceptional lagrangian singularities so that (A; R) is stable under small perturba- 
tions. Let also S > 0. Then there is a sequence of tori A-^'^^ , for J' = J'^{h) satisfying 
\k\h < , such that H has quasi-modes {uk{x, h), Ek{h)), \Ei^{h) — E\ — 0{h^), in the sense 
{H{x, hDx) — Ek{h))uk{x, h) — 0{h'^). Morover, these quasi-modes are of maximal density, 
in a /^'^-neighborhood of A. 

Remark 2.5: In the case H contains a subprincipal symbol and equation (b.29) holds, the sub- 
principal form makes sense on A, only when it is invariant by Xh- But again, the diophantine 
condition and our special action-angle variables allow to extend these constructions to the 
case where Xh is no longer tangent to A, provided we solve (b.29) as in (b.31) instead of 
(b.25), and modify Maslov quantization condition as in (c.4). 

3. Jacobi-Maupertuis correspondance. 

We will apply the constructions above to Hamiltonian H{p, x) = + V{x) near energy 
E — 0, where V{x) = —S/g{x). The "reference Hamiltonian" Ti. = g{x)p'^ denotes an 
(covariant) metric on T*M, for which A C {7i = S} is an integral lagrangian torus with 
irrationnal frequencies uj. The Hamilton vector fields are related by Xj-i = g{x)XH, so that 
if T parametrizes the integral curves of X-j-i and t (the new time) those of Xh, we have the 
relation dt = g{X'^{T))dT. This is Maupertuis-Jacobi correspondance. 

a) Construction of quasi-modes. 

First we show that, assuming a diophantine condition on u, the flow on A for H is 
again quasi-periodic. Expand g e C°° as a Fourier series on A R'^/27rZ'^, g(^X^{(p)) = 
{g) + ^ c^e*"*^, Cn rapidly decreasing, we have : 

n6Z''\0 

Proposition 3.1: Assume a diophantine condition on a;. Then Maupertuis-Jacobi corre- 
spondance induces a reparametrization of A by 

^ : (fio ^ (fi = ^{(fio) = (fio+ (^fifo) 

where / is a smooth, scalar periodic function with (/) = 0. The motion on A induced by this 
reparametrization is quasi-periodic with frequency vector a; — uj/{g), and (for simplicity we 
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denoted g o X° by g) : 

9$ (a) 

(3.1) det|-M = ^ 

ofo 9 

Proof. Integrating dt = giujr + (po)dT we get 

t^tQ + {g)T + FiujT + ^o) - F{ipQ) 

where F is obtained by integrating g term by term, and to is an arbitrary constant. We 
have (F) = 0, and F is again smooth due to the diophantine condition. Choose to = -^(v^o), 
and change variable z = r — t/{g). With the new frequency ui — u}/{g), the latter equation 
rewrites as 

(3.2) h{z,i(;t) = z + F{il;t + ^z)/{g) = 

where we have set ifjt = f^o + ojt. Equation (3.1) holds when z = zq = —to/{g), ifjt — 
V'to = Vo + '^to, and moreover ^ = g{'i(jt + u)z)/{g) > 0. So implicit function theorem 
shows that the solution to (3.2) is of the form z = f{ipt)-, i^t G A, where / is periodic and 
(/) = 0. So far we have shown that Maupertuis-Jacobi corrcspondance can be written as 
ipo + ujT = (^0 + ujt + ujf{(^o + ivt). To get the Proposition, it suffices to observe that the 
trajectory t ^ <fo + cot is everywhere dense in A, since u is again irrationnal, and moreover 
^((^ + ujf{ip) = I + ujf'{i-p) is of maximal rank, because 

det(/ + cJ/'(^)) = M 
9 

which brings the proof to an end. ^ 

Note that uj and u) are simultanously diophantine. We can now construct quasi-modes 
for H supported near A using the procedure set up in Sect. 2. On A, Maupertuis-Jacobi 
corrcspondance will only induce by (3.2) a change in the densities, namely 

(3.3) det —— = ^ det - — , det —— = ^ det - — 

dcp g dipo dip g dcpo 

and also in the mixed representation. 

On neighboring tori A"^''^ , we have (J', ip') — k'^{J^ ip) = k.]^^{J, ^{(po)), so that we can 
still consider (p', for given J, as a function of (po alone. Let us summarize our result so far in 
the : 

Theorem 3.2: Let Ti. = g{x)p^ be a (not necessarily integrable) covariant metric on T*M 
have an invariant torus A with T-C\a = ^, and / be the actions as in (1.1). Assume that 
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on A the Hamiltonian flow of H. has diophantine frequencies a;, and morerover A is free of 
exceptional lagrangian singularities so that (A; R) is stable under small perturbations. 

On the other hand let H{p^ x) = + V{x), where V{x) = —E/g{x), be the hamiltonian 
on T*M given by Maupertuis-Jacobi correspondaance. Let also 5 > 0. Then there is a 
sequence of tori A-^', for J' = J'^{h) satisfying \k\h < , such that H near energy E = 
has quasi-modes {uk{x, h), Ek{h)),\Ek{h) — E\ = 0{h^), supported on A-^'^^, in the sense 
{H{x,hDj;) — Ek{h))uk{x,h) = 0{h?). Morover, these quasi-modes are of maximal density, 
in a /i^-neighborhood of A. The quasi-eigenfunctions u{x,i'^h) = Uk{x,h) can be expressed 
in term of "Maslov canonical operator", with the help of formulas (2.15), (2.19), (2.20). In 
particular where A is projectable, we have : 
(3.4) 

where (po G T''. So the half-density | det ^^g'^/'^ ^ 1""*"^^? for given J', can be written as a 
function of (fio alone. 

Instead of H{x,hD^) = -A'^ + V{x) we may consider a h-PDO with symbol H{p, x; h) = 
+ V{x) + hHi{p,x) + ■ • ■. The introduction of a sub- principal symbol can be handled 
similarly, due to Remark 2.5. It modifies | det ^^g^/'*^ ^ by a factor 

eM-iGi'fi%j',^')=n-\jM^o)) 

as in (b.31), and Maslov index by the vector (crH)|AJ',Jv(l, • • • , 1) as in (c.4). 

More generally, this procedure may be applied whenever Maupertuis-Jacobi correspon- 
dance relates the flows of two classical Hamiltonians H and H on a lagrangian manifold A, 
H not necessarily of the form p'^ — l/g{x), as is the case for the shallow water waves disper- 
sion relation of Sect. 4. d. Actually, all what we need is to construct an oscillating principal 
symbol for lagrangian distributions microlocalized on the family A"^ and this can be done 
following the standard procedure we recall in Appendix b. 

b) Remarks on isospectrality. 

There naturally arises the question of comparing the quasi-modes microlocalized near 
A of the Schrodinger operator H = — /i^A — E/g{x), as given in Theorem 3.2, with those 
of the Laplacian l-L{x,hDx)hDxg{x)hDx on M, in other words : does Jacobi-Maupertuis 
correspondance imply (partial) isospectrality ? Theorem 3.2 gives all quasi-modes of H 
supported in a microlocal neighborhood of an invariant torus A for Ti. Here we address the 
converse problem, i.e. how to construct quasi-modes for the Laplace operator associated with 
H. (that we may write in a symmetric form as 7i(a;, hD) — {hDx, g{x)hDx) see e.g. [Mat], ) 
with quasi-eigenvalues close to the Ek{h) ? 
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We shall assume here that Ti. = g{x)p^ is integrable, so that in global action-angle 
variables, Ti. = H{I) = {uj{I),I), where / are as in (1.1), and uj{I) = Choose I = Iq 
as a base point for the action-angle coordinates (J', (f') we have constructed so far for i?, 
and consider 'H{J') = {oj{J'), J'). Under some generic non degeneracy condition, the integral 
manifolds oi H{J') are close to those of 'H{I), for small enough t'. By definition, the action 
integrals on the lagrangian tori A-^' (integral manifolds for H) verify (1.1), and the A-^' are 
close enough to the A-^ '^ (integral manifolds for H), although in general A-^' ^ A-^ '^ if 
t' ^ 0. It is still possible to compare the action J pdx along A with the action J pdx along 
A-^' as in Propositions 2.1-3. We can also extend the WKB constructions leading to Theorem 
2.4 in the case of the Schrodinger operator H to the Laplace operator H{x, hD), and the fact 
that X-j-i is tangent to A"^' makes them easier. To get Maslov quantization condition on A-^' 
it suffices to choose J' = J'^{h) (instead of J' = J'k{h)). Since J'^{h) = J'k{h){l + 0{h)), this 
gives easily the required (approximate) isospectrality result. 

4. Some applications 

We start by examining the simple case of a conformal, radially symmetric metric. Then 
both geodesic flow and hamiltonian flow (with potential) are integrable, and the constructions 
above are not truly necessary, but this can help to understand more complicated situations, 
which we consider next. Finally we study the Hamiltonian associated with the dispersion 
relation that arises in the shallow water waves theory. 

a) Radial potential and related Zeeman effect. 

The metric in cotangent space T*R^ we consider is given by 7Y(p, x) = p'^g{\x\). In polar 
coordinates x = (r, 6) with dual variables p — {p,ip), the kinetic energy reads p^ — + ^ 
and Hamilton equations take the form 

f = 2pg{rl p= V(r)(p^+^)+2^V^ 

(4.1) ^ ^ 

(Because of the radial symmetry, this is equivalent to consider the Hamiltonian with poten- 
tial.) This example is of course well known (see [Ze], [HiSjVu] and references therein) but we 
want here to give an explicit parametrization for (r, 9) in action-angle variables, from which 
we deduce the configuration of caustics. We assume the energy surface Eg 

(4.2) g(r){p^ + ^)=E 

to be compact, and look for the lagrangian foliation of Eg. The last equation (4.1) shows that 
ip = ipo is a constant of the motion, so the moment map is of the form m = {g{r) (p^ -|- ^) , ip). 
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Let us find the actions variables and the fundamental frequencies u = {u!i,u!2) on the torus 
A = {m = {£, ipo)}- We have 

(4.3) (P, dX) = pdr + ijde 

Assume there is ro > such that ipQ^^^ = £. Then C = {r = ro} is a caustic of A 
(turning surface) and for r = tq, p = 0. For a general Morse function g{r), there are 
several turning surfaces, delimiting annuli of radius ri < • • • < in the a;-space. Choose 
an annulus D = {ri < r < r2} in the classically allowed region, and call A the connected 
component of the integral manifold lying above D. As a first fundamental cycle 71 in A, we 
take r = ro, (ro = ri or r2), V' = V'o- Integrating (4.3) along 71, we find the first action 
variable Ii = I^yj^iP^dx) = V'o (angular momentum). Let 72 the second fundamental cycle 
be given by ^ = ^O; '^i ^ r < r2, this yields Clairaut integral 

Substituting t/^q = /i, since the turning points ri, r2 depend also on Ii, we find that I2 is of the 
form I2 = f{Ii) where / is a smooth function (recall that £ is held fixed.) The hamiltonian 
reduces to H{h, f{h)) = £, and differentiating this relation we get the rotation (or winding) 
number 

Next we notice that the 2 first equations (4.1) relative to the (r, p)-projection of motion are 
decoupled from the other equations for the ^-projection. So integrating these and using that 
A is a torus of periods 2tt, we find that 

(4.4) r = r{u2t + cp^2), P = p{uJ2t + ^l) 

for some constant ip2- Now we use the third equation (4.1) which can be integrated as 
6{t) = 2il)Q J* ^^\r=r{L02s+^%)ds ■ We expand as a Fourier series in a;2-s + V'25 integrate 
term by term, which gives 

(4.5) e = 6{L02t + ^l)+uJit + ^p\ 

Here (/?° is a constant, and = 2'0o(-^^), the 0:th order Fourier coefficient, ((■) denotes the 
time average over [0,27r].) 

We look now at the caustic and determine the focal points. Note that because of con- 
servation of momentum, to any x inside the annulus, correspond exactly 2 values of p. We 
have 

ax ^ /^cos^-rsin^^ -rsin6'\ ^^^^/^cosV 
d(p y^sine + rcose^ rcoso)' dip y^smip 
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From (4.1) and the relation ^ = <^ig^ + <^2a^ on the torus A, we recover the fact that 
the focal points are given by = 0, i.e. p = 0. Moreover, the condition that (C, B) is 
everywhere of rank 2 (see Leemma a.l ) yields ^ when p = 0. The caustic r = tq 
consists here only in points of Tq^ and the degeneracy is due to the radial symmetry of the 
metric. A particular case, with zero angular momentum, was studied in [0-de-AlHa]. 

Remark: It is known that due to the positive homogeneity of Ti with respect to the action 
variables, Euler formula gives cui/i +UJ2I2 = ^ (see [Ze]). Let us give another simple proof of 
this. Considering the Lagrangian formulation associated with 7Y, we see that the energy of 
the (unit mass) particle moving along the geodesic flow is 
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where ^ denotes limT-»oo Jq ■ Substituting x = {(jj,d^)x, we get £ = ^^oJj 4 {-^^,p)dt. 

Now if the frequencies are rationnally independent we can replace the time average by the 
average over angles variables, so that by the definition of Ij : 

^ = Zl'^i^^Tr)-^ / {-^x,p)difi - ^ujjlj 

Consider now Maupertuis- Jacobi correspondance allowing to switch from Ti. to the hamil- 
tonian H = p^ — As a semiclassical example, we consider the Schrodinger operator with 
potential —u{r) : 

H - - — —r— u(r) 
r dr dr d9^ 

We separate variables expanding as a Fourier series in 9, and rescale by ^/r, which leads 
to a family of 1-d semi-classical channel operators of the form Hn = —h?-§:p2 + Vnir) with 
effective potential Vn{r) = — u{r). Thus the angular momenta are given hy ip = 

ipn = h{n'^ + 1/4)-*^/^, and the turning surfaces r = tq at energy E by r'^{E + u{r)) = i/jq. So 
because of the centrifugal barrier, the energy shell for can have a bounded component 
< ri < r < r2 in configuration space, and another unbounded component r > r^ (consider 
e.g. u{r) = 1/r + r > and E = 0. ) A torus A corresponding to the angular momentum t/j^ 
in the compact component of the energy surface, carries a quasi-mode for suitable E = E{h), 
and E{h) turns in general into a resonance for hamiltonian H, whose imaginary part has to 
be computed. 

A famous example is given by the dynamics in R"^ of an electron in a Coulomb potential 
— 1/r , and interacting with a small, constant magnetic field of strenght B (the constructions 
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above should be extended to the 3-D case). In suitable units, and cylindrical coordinates 
(p, 9, z), with dual variables ipp,Pe,Pz) the Hamiltonian takes the form : 

H = {pI+pI + ^) + Bp' - ip' + z')-'/' 

The Coulomb singularity can be removed by introducing Jacobi or "parabolic" coordinates, 
see e.g. [Ts] . The Rayleigh-Schrodinger series for resonances near E = 1/4 is Borel summable 
and can be computed perturbatively around B = (Zeeman effect). This (non semiclassical) 
problem has been completely solved in [HeSj] , and lies far beyond the scope of our analysis. 
Nevertheless, we would like to point out that resonances arise from connecting real tori (one 
of them is actually a cylinder, and lies over an annulus of the form Q = {ra < p < +oo}, ) 
through tunnel cycles in the classically forbidden region, see [DoSh] for a related problem. 

b) The Kerr-Newman space time. 

We give another example of a metric for which we can find some invariant (degener- 
ate) tori and reparametrize the natural phase space variables with action-angle variables by 
the procedure above. Namely, consider the geodesic fiow associated with a simplified Kerr- 
Newman metric (see e.g. [Ch], or [Ha] and references therein, ) for which the black- hole 
angular momentum per unit mass is a = (Reissner- Nordstrom space-time), M its mass and 
Q its charge. The Hamiltonian reads : 

1 4 2 *2 

2r2^ A sin^^ ^ 

The variables (t, r, r, p, 6', 6'*, (/?*) G T*R^+^ are denoted collectively by {P,X), and A = 
A(r) = -2Mr + (5^ is a second order polynomial with constant coefficients having 2 real 
and positive roots ri < r2- The region ri < r < r2 is called Boyer-Lindquist block Bu. The 
equations of motion are of the form 

; r2 . ^ . A . dn 
/A r,\ A or 

2H sm 6* sm 6* 

It is known that there are 4 integrals of motion, given by 

t2 4 2 

p = 2n,E = T,L = -<f*,)C = e*' + = — — - Ap^ - pr^ 



sm 



9 A 



The energy surface is of course non bounded, but has some bounded connected components 
in the space-like region t — to,T — 0, we can call "unphysical" . The moment map restricted 
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to T = is of the form m = (p, L,/C). Let us find the actions variables and the fundamental 
frequencies uj = (a;i,a;2,(^3) on the corresponding torus A. We have (P, dX) = rdt + pdr + 
9*d9 + (p*dip. 

As a first fundamental cycle 71 in A, we take r = tq & {ri,r2}), (fi* = —L. Integrating 

* 

(4.7) along 71, we find that dcp and d9 are related by dcp = QT^^dd so the first action 
variable takes the form : 

= -1 / {P,dx) = ^ r(}c - J^y'/'de 

which makes sense only if L = 0. So the dimension of A is 2, and Ii = VK.. 

Let 72 the second fundamental cycle be given by ^ = '^i ^ r < r2, this yields another 
Clairaut type integral 

where we have made use of £^ = 0. 

The last components of the degenerate (3+l)-torus int = to is then 3 = ^ Jq^ (p*d(p = 
0, and ujs is undefined. It is then clear from the Hamilton equations (4.7) that the description 
above of the motion for the geodesic flow with a radial metric extends to this case without 
change. In particular, formulae corresponding to (4.4) and (4.5) hold. We shall not investigate 
further the case of degenerate tori (see [BeDoMa] for a general discussion.) Other "solvable 
models" are considered in [DuFoNo,Chap.7]. 

c) Liouville metrics. 

These metric have been extensively studied ( [Fo] , [BNe] , [Se] , [Mat] , . . . ) for their nice 
property of being square integrable, in the sense that the geodesic flow has another quadratic 
integral, i.e. a quadratic polynomial in p. 

We call Liouville metric on a domain Q C a conformal metric given in T*Q by the 
hamiltonian T-[{x,p) — g{x)p^ with a quadratic (in momentum) second integral. 

We call model Liouville metric on VL the Liouville metric where 

(4.8) g{x) = {u{xi) + v{x2)) 

It is straightforward to check that, besides 'H{x , p) — E , oxv additional integral of the motion 
is given by 

(4.9) fip^x)- 



u{xi) + V{X2) 
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A classical result of Birkhoff (see [Fo,p.291]) says that in certain appropriate coordi- 
nates any 2-D Riemannian metric whose geodesic flow has a additional integral quadratic in 
momenta is (locally) a model Liouville metric. 

An interesting case is when g is globally of the form (4.8) on Q,. Thus the hamiltonian 
system decouples completely into two 1-d hamiltonians with a potential pf = u{xi)E-\-F^p2 = 
v{x2)E — F, and moment map m = {E,F). The corresponding semiclassical Schrodinger 
operator is also completely integrable, in the sense that the spectral problem reduces in 
finding the joint spectrum of 2 commuting 1-D Schrodinger operators, see e.g. [KMS] for a 
study of the spectra of Laplace-Beltrami operators on a Liouville surface, and also [DoSh] for 
tunneling effects. 

Moreover, square integrable conformal metrics when Q is an annulus, are globally, con- 
formally equivalent to model Liouville metrics in the following sense : 

Proposition 4.1 [Mat]: Let 7i(x,p) = g{x)p'^ be a smooth conformal square integrable 
metric on a compact domain O of the x-plane diffeomorphic to an annulus, and A C T*0 be 
an integral lagrangian manifold for the hamiltonian flow such that tt : A — > O is surjective 
(A has no "holes" , i.e. there is no forbidden region above ) Then there exists a bijective 
conformal map a : O — * O (depending on A) where O is an annulus of the x-plane, and such 
that the image of Ti through a becomes a model Liouville metric. 

One can also compute action integrals / and corresponding frequencies u as in Part a). 
In particular we have the relation ujili + UJ2I2 = 2i?. At the level of operators, we use Weyl 
/i-quantization of the form T-C{x,hD) = hDxg{x)hDx, so that the diffeomorphism a takes 
H{x,hD), denoting by y the new space coordinates, to 7i{y,hD) = hD\Wa\^g o a{y)hD. 
This operator is of special interest in the theory of surface water waves captured by coasts 
and bottom relief irregularities. 

Consider now the following semiclassical "Maupertuis problem" : let H{x,p) = g{x)p^ 
and A C T*Q, be as in Proposition 4.1, with 7Y|a = 1. Let H = -h'^A-l/g{x), and A = A{x) 
be the unitary Fourier Integral Operator (FIO) on L'^{Q) implementing a (suitably extended 
to the whole space). Then H = A*HA is a self-adjoint h-PBO on L^(l], iVapcZy) with 
principal symbol Ho^y^p) = p^ — l/g{y) (here g = g o a is the model Liouville metric on Q 
given in Proposition 4.1. ) There is no sub-principal symbol. Moreover, A is an invariant 
lagrangian manifold for X~. So we may apply Theorem 3.2 to construct quasi-modes for H. 
In the following subsection we shall introduce another Maupertuis correspondance for this 
geodesic flow. 

d) The shallow water waves. 

Let Q C (the "basin") be the complement of a bounded set (the "island"), so that 
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Q is difFeomorphic to an (unbounded) annulus. In the shallow water waves theory [DoZh] 
one constructs a /i-PDO with symbol L{p, x; h) = Lq{p, x) + hLi{p, x) + ■ ■ ■. The principal 
symbol of L (dispersion relation) on T*Q reads : 

(4.10) Loip,x) = |p|(l + /xp2)tanh(D(a;)|p|) 

where D is a smooth positive function, describing the depth of the basin, and fx > a, small 
parameter. Quasi-modes of L represent the waves trapped by the island. By a Maupertuis- 
Jacobi correspondance, we shall construct trapped modes for L at energy E from lagrangian 
tori of Liouville metrics as above, by choosing the depth profile in an appropriate way. Since 
Liouville metrics are actually generic among metrics with a quadratic second integral [BeN] , 
this result could be of practical interest. 

For simplicity we shall assume /i = 0, but it is not hard to see that our constructions 
may be carried out perturbatively, due to implicit functions theorem, for > small enough. 
Note that, if D{x) is positively homogeneous of degree 1, then integrability of Lq{p,x) = 
\p\ tanh(D(a:) |p|) at energy 1 implies integrability at all energies E > (make the symplectic 
change of coordinates {x',p') = {Ex, E~^p). ) In general case however, Lq{x,p) has no reason 
for being integrable at any energy E, but finding suitable profiles that give raise to invariant 
tori can be achieved by the following Maupertuis correspondance. 

Namely, multiply relation Lo{p,x) = E (for ;U = 0) by D{x), and set z = \p\D{x), 
(j){z^) = ^tanh^;. We denote by ip : R+ — > R+ the inverse of (j){z^), which is a smooth 
increasing function, so that (4.10) can be rewritten as 

(4.11) n{x,p) = g{x,E)p' = l, g{x,E)- ^'^^^ 



i;{ED{x)) 
We have the elementary : 

Lemma 4.2: On the energy surface H{x,p) = 1, the Hamilton vector fields Xlq and X-h 
are related by X-^ = di^^^^/^ed{x)) -^h, and this induces a reparametrization of the integral 
curves of Xjj{X{t), P{t)) given by Maupertuis- Jacobi correspondance 

^ ^ D{x)ip'{ED{x)) 

So we need look for depth profiles D{x) (depending also on E, ) for which g{x,E)p^ 
is integrable. Clearly the function G : R+ [0, 1[ defined by G{D) = ^^j-,^ is a smooth 
increasing function, satisfying G{D) = D{1 — + ■ ■ ■) as D ^ 0, and lim^^oo G!{D) — 1, 
so the depth profile D{x) = D{x, E) can be chosen as the solution of g{x, E) — g{x) = 
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E '^G{ED{x, E)), where g{x) is one of the square integrable metrics given in Proposition 
4.1. 

Now if the geodesic flow H{x,p) = g{x, E)p^ is integrable and admits invariant tori, the 
procedure set up in Sect. 3 can be extended to construct quasi-modes for Weyl /i-quantization 
L'^{x, hD) of hamiltonian L{x,p; h), for which eigenfunctions, or quasi-modes, describe again 
trapped waves near energy E. Note however that L is a h-PT>0 of more general type than 
a Schrodinger operator, but all our constructions, and in particular the results of Appendix, 
apply to this setting, as we pointed out after Theorem 3.2. Thus we can state our main result 
as : 

Theorem 4.3: Consider the h-PDO with symbol L{p,x;h) — Lq{p,x) + hLi{p,x) + ••• 
with Lq as in (4.10) defined on an open set Q, C diffeomorphic to an annulus. For all 
£■ > 0, and all smooth Liouville metric T-t{x,p) — g{x,E)p^ on T*Q,, there exists a "depth 
profile" D{x) ~ D{x,E) on fi, such that the pair of Hamiltonians (Lq^H) have the same 
trajectories at energies {E, 1), both parametrizations being related by (4.12). In particular, 
T-C{x,p) being integrable, {Lo,7i) have a common family of lagrangian tori A for energies 
{E, 1), and Theorem 3.2 allows to construct a sequence of quasi-modes for L{x, hD^, h) near 
energy E. 

Appendix : Regular WKB constructions and Maslov theory 

In this Section we recall some well known facts about semi-classical quantization by 
standard Maslov theory ; we refer to [Arl,2], [So], [CdVl], [Hol,2], [Laz], [Le], [Li], [Dul,2], 
[GuSt], [Iv], [We], [MeSj], [BaWe], [DoZh], [deG]...for more details. In Part A we give a 
rather superficial account for singularities of a lagrangian manifold diffeomorphic to a torus, 
and paramerized globally by angle variables e T'^ (essentially when d = 2) ; then we 
derive a simple formula for Maslov index. In Part B we recall some known facts about WKB 
theory, when the tangent space to the lagrangian manifold contains Xh-, as a preliminary 
step to Sect. 3 for the perturbative case. In Part C we combine these results to get Bohr- 
Sommerefeld-Maslov quantization condition. 

A. Maslov index and singularities for a lagrangian immersion of the torus. 

Let M be a smooth Riemannian manifold, (we shall think essentially of M = R*^ or the 
flat torus M = T'^ in the case of the geodesic flow on a Liouville surface, ) and l : T*^ — T*M 
a smooth lagrangian immersion, we parametrize A = t(T^) hy p = P{(f),x = X{(f). To 
comply with usual notations, we identify already T'^ with its image through l, so that the 
Lagrangian immersion reads now l : A T*M. 

We begin with a short description of lagrangian immersions. Contrary to local la- 
grangian singularities, the global geometry of lagrangian tori, especially those arising from 
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integrable systems is not so commonly discussed (see however [0-de-AlHa], [DoZh] and ref- 
erences therein. ) We call focal point a point y E A where tt* : TA TM is singular, 
and caustics the projection C of the set of focal points onto M. We introduce the jacobian 
matrices B = ^ and C = and recall the following standard result (see e.g. [DoZh]) 

Lemma a.l: With the notations above, we have : 

1) the matrix {B, C) is of rank d. 

2) the matrices ^CB, BC~^ and CB~^ are symmetric (whenever well defined, otherwise 
replace e.g. B~^ by dei{B)B~^ . ) 

3) C ±iB is non degenerate. 

Proof. The 2 first statements easily follow from the definition of a Lagrangian manifold, and 
the 3:rd one from computing the hermitian scalar product (^|(*5C — *CS)^) when ^ belongs 
to the kernel oi C ±iB. X 

Taking variations, we see that tangent vectors {5X, 5P) e TA are given by 5P = B{y)5(f, 
6X = C{y)5^, so 6P = BC-^{y)5X when TyA is transverse to the fiber, or 6X = CB-^{y)5P 
when TyA is transverse to the zero section of T*M. So (f)"{x) = BC~^{y), when A is 
parametrized by the 1-form dcj) = L*{pdx), while il'"{p) = CB~^{y), when A is parametrized 
by the 1-form dip = i*{—xdp). 

For simplicity restrict now to the 2-d case, and let y be a focal point. By property 1) of 
the Lemma, either B{y) is of maximal rank, or both C{y) and B{ii) are of rank 1. 

In the first case, either C{y) is of rank 1 (which, under some transversality condition, is 
the generic case, and corresponds to the simple fold A2, ) and A near y intersects (tangentially) 
the fiber VM (the vertical space of T*M) along a line, so C is of codimension 1, or A is tangent 
to the fiber of T*M at j/, and we call y an umbilic. The set of umbilics denoted by is 
generically discrete, and thus does not contribute to Maslov index. 

In the second case there arc A, /i G R \ such that -f^ = A-^^ and 4^ = u-i^ at y. 
Identifying the off-diagonal terms of *i3C, which is symmetric by property 2) of the Lemma, 
we find that cither X = jj or ± (in the sense of the cuclidcan structure on T*M, since 
M is flat. ) But A 7^ /U since otherwise the complex matrices C ± iB would be degenerate, 
which violates property 3) of the Lemma. So and are colinear, and orthogonal to 
both 1^ and 1^. 

We call y a cusp. The corresponding lagrangian singularity is generically the Whitney 
crease A3, and we denote the set of cusps by J-'q. 

Assume now A C E^; where is a level set for an Hamiltonian of the form H = 
p^ -\- V{x). Differentiating, we get 2(P, + (VF, ^-) = on A at y, j = 1,2, and since 

II, these relations imply |g- ± VV{X) and P ± j = 1, 2. So if both C{y) and B{y) 
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are of rank 1, the vectors P, and are colinear, and all orthogonal to S/V, and 

dP 

Generically (i.e. without additionnal symmetries) the set of such focal points y is 
discrete, or even empty, and like umbilics, do not contribute to Maslov index. The same 
argument applies \i E. = p^g{x) (geodesic flow) or H = {pi — x^)^ + (p2 + x\)^ + Vix). See 
however Sect. 4 for the case of rotational symmetry. 

Example a.l: The harmonic oscillator H = p"^ + Vo{x), Vo{x) = u^xl + lo\x\. Although 
standard action-angle variables are singular at energy i? = 0, they are often best suited for 
describing the singularities at caustics. For all £^ > 0, there are 4 (hyperbolic) umbilics, 
located at the vertices of the caustics (a rectangle). For small energies £' > 0, these umbilics 
are stable under perturbations of the form V{x) = Vq{x) + C(|a;|"^). 

Example a. 2: The magnetic hamiltonian in of the form Ha = {p—A{x))'^+V{x) (written in 
the minimal coupling setting) with V as above, and A{x) = 0{\x\'^) near x = 0. Then umbilics 
(when A{x) = 0) metamorphose into cusps, and singularities are precisely Whitney creases, 
see [KaRo]. For more general magnetic hamiltonians in T*M, M = R"^, it is convenient 
to change the minimal coupling formulation, and introduce the non-degenerate (magnetic) 
2-form as — o — B{x) with B{x) = dA(x), which makes of T*M a new symplectic manifold. 
When V{x) = ^^Vj(xj), for smooth potentials Vj{xj), it is easy to see that Ij — {pj — 

3 

Aj{x))'^ + Vj{xj) Poisson commute for ctb, and thus can be chosen as action- variables. 

Example a. 3: A radially symmetric metric Ti = g{\x\)p'^ on an annulus < ri < r < r2- In 
general, the whole caustics lies in J^o (see again Sect. 4). (This emphasizes the peculiarity of 
the harmonic oscillator among integrable systems. ) 

Example a. 4- The geodesic flow on an ellipsoid is conformally equivalent, in the upper-half 
space, to a Liouville metric H = g{x)p^ on a rectangle =]0, Ti [x]0, T2[. There are 4 umbilics 
going to the vertices of i?, and the induced metric near an umbilic takes the form, in local 
coordinates l/g{x) = G{x\) - G{-x1) = (a;f + xl)A{xl - x\, 2x1X2), where A{0, 0) > (see 
[CdVVu,Sect.3] and references therein. ) 

Now we can pass to Maslov theory and consider first the linear situation. Recall the real 
phase space T*R^, in addition with its symplectic and euclidean structures, is endowed with 
the complex structure J{p,x) = {—x,p), ; we identify T*IV^ with and its zero section 
Rp © with the totally real subspace R^. 

Let C{2d) be the lagrangian grassmannian of T*R'^, we know that U{d) (the unitary 
group on C^) preserves all these structures and acts transitively on C{2d). We can identify 
C{2d) with U{d)/0{d), 0{d) (the orthogonal group on R*^) being the stabilizer of a given 
R'^ fti A e L{2d) under the action of U{d). Equivalently, if A = Wi(R'') = W2(R'^) (with 
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the identification above of the zero section, ) Uj G U{d)^ then U2 =U\0^ O & 0{d). Since 
det^C = 1, det^Wi = det^W2 only depends on A. The map det^ : U{d) — > induces 
a fibration C{2d) — > S-^, with 1-connected fiber SU{d)/SO{d), giving an isomorphism of 
fundamental groups : 7ri{jC{2d)) 7ri(S-'^) Z. Passing to homology and dualizing, we obtain 
a natural homomorphism : H^{S^; Z) H^{C{2d) \ Ti). The image ji = jid oi the canonical 
generator of H^{^^-^ Ti) under this map is called the universal Maslov class. Explicitely, along 
the path C{2d),t ^ At = Wt(R^), Ut e U{d), we have 

I 

(a. 3) diJ,{At) = d/it = — d(argdetZ^^^) 

ZTT 

Consider now the symplectic vector bundle tt : T*M — > M over a manifold M, with a 
lagrangian subbundle A — > T*M. In general, the automorphism group of T*M does not act 
transitively on the lagrangian subbundles, but a pair of transverse lagrangian subbundles can 
be related as follows. Let <, : A — > T*M be a lagrangian immersion, then the symplectic 
vector bundle l*T{T*M) over A has 2 lagrangian subbundles L = i*TA and L2 = l*VM (the 
"vertical" subbundle). The same consideration holds when replacing L2 by the zero-section Li 
of T*M. Maslov class a = ij,a,i. oi (A, t) is defined as /i^^a = l^{L, L2) e -^''^(A; Z), which is the 
pull-back of the universal Maslov class ji by Gauss map G : A — > C{2d)^ y 1— > G{y) = i^TyK. 

For y e A, let U{y) be the unitary transform generating the lagrangian plane TyK. 
Possibly after changing djit into —djit in (a.3) it is clear that, at least for d = 1, U'^{y) is 
given globally by Cayley transformation 

(a.4) U\y) = {C + iB){C-iB)-\y), d=l 

Example a. 5: In case A is the product of d circles + = Ej parametrized by (pj^xj) = 
{y^Ej siinpj, cos (fij), (a.3) and (a.4) give by integration Maslov index a — /x(L, L2) = 
(2,---,2). 

We still show below how (a.4) holds "at the infinitesimal level" and allows to compute 
the jump of Maslov index (u(L, L2) G H^{A, R) passing a caustics. See also [Ma2], [DoZh] for 
another and more explicit derivation of /u(L, L2). 

Moreover, from the discussion after Lemma a.l, the RHS of (a.4) equals {I + iBC~^){I — 
iBC~^)~^{y) when A is transverse at y to L2, and —{I — iCB~^){I + iCB~^)~^{y) when A is 
transverse at y to Li. Thus, because BC~^{y) and CB~^{y) are symmetric, except perhaps 
on an exceptional set which we identified with JFq in the 2-d case, the RHS of (a.4) is well 
defined as a unitary matrix acting on lagrangian planes, i.e. the fibers of TA. So we are in the 
situation of [Ho2,p.l56], up to the sign in front of — (/ — iCB~^)(I + iCB~^)~^{y), which we 
may ignore if d is even (since only det U^{y) matters), or if we forget about orientation. So we 
writeU^iy) = {I + iBC-^){I - iBC'^y^y) (resp. U^{y) = {I - iCB-^){I + iCB-^)-\y). 
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Following [So] (see also [GuSt]) we shall identify these unitary matrices with the la- 
grangian planes on which they act, and introduce also the covering space TA of TA by 
adding a new (angular) variable 6 which stands for arg det U. So in TA we have the freedom 
to consider a path of "fictitious" lagrangean planes {Lt)o<t<i which are not tangent to A, 
except for t = 0,1. A similar argument, consisting in interpolating between C + ieB and 
C — ie B, was used in [DoZh]. Now we argue as in [Ho2], and to start with, restrict to the 
case where Cy^ has codimension 1 in A and B{y) is of maximal rank d near j/q- 

Consider first a path (7i(t), 7i(t)) G TA, — eo < t < cq, with 71 (t) G T^^(4)A = T{t), and 
such that 7i(0) belongs to the fiber. So in the terminology of [GuSt,p.l20], t(— eo) and T(eo), 
lie in different components of the projective lines jCl^ (corresponding to the fiber L2) and 
jCl2 (corresponding to the zero-section Li). We lift the path (71, 71)) in TA as (71, 7i). Then 
along this path, 5(y)|j^=^j(i) is always invertible, and W^(j/) = {I — iCB~^){I + iCB~^)~^{y). 
The square of the determinant of Ut is Ylji^ ~ + where Pj{t) are the real 

eigenvalues of CB~^. Since arg(l — i(3j{t)) = — arg(l + if3j{t)) (we choose argz so that it 

vanishes for 2; > 0) we get d/it = ^ ^ ai'g(l — iPji^o)) — arg(l — iPj{—eo)). 

j 

Choose next a path (72(^)5 72 (^)) ^ with 9 between argdetW(7(eo)) and arg det 
W(7(— eo)), and whose projection starts at r(eo) and ends at r(— eo). We use instead V(f{y) — 
{I+iBC-^){I-iBC-^)-\y). The square of the determinant of Ut is now Yl^{l + i/ Pj{t)){l- 

i/Pj{t))-\ and/~//t = i5;]arg(l+V/3,(l))-arg(l-V/3,(0)). Since for real /3 ^ 0, 

j 

arg(l - ip) - arg(l + i/P)) = -7r/2 if /? > 0, = 7r/2 if /? < 0, we find : 

cfdfit={f + hdfit^l{sgn{BC-'){eo)-sgn{BC-'){-eo)) 

Here f denotes the integral over the loop 7 = 71 V 72 in TA, and we have used that for t 0, 
BC~^ and CB~^ are well defined and sgn5C~^ = — sgnC5~^. This gives the required 
formula 

(a.6) (7,«) = ^(sgn(SC-i)(eo) -sgn(SC-i)(-eo)) 

When J-'o is a discrete set it doesn't contribute to Maslov index, so formula (a.6) gives 
generically the jump of Maslov index passing a caustic. But in some cases as the radial 
symmetric potentials in d = 2 considered in Sect. 4, there are caustic sets where both B{y) 
and C{y) have rank 1, and the constructions above fail. The argument of [DoZh] covers most 
of these exceptional cases, but we can still adapt this one as follows. 

Actually, since we can construct a subbundle transversal to both L2 = l*VM and L = 
t*TA, we know that near any yo = (po, xq) e A, there is a local diffeomorphism x (in the 
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variable x alone), such that A is parametrized by a non-degenerate phase function of the type 
S{p,x) = (p,x) + H{p), i.e. A = {(p, a;) : dpS{p,x) = 0,p e neigh(po)}- Moreover, H"{p) is 
symmetric and invertible. 

Example a. 6: Let A C T*R be the circle y'^+r]'^ = 1, after the symplectic change of coordinates 
y = x(x) = — 1,?7 = *(x'(^))~^C — near (2/,?7) = (0, 1), A rewrites as a; = ^ — 1 + 

When A is of the form {P{(p), X ((/?)), </? G T^, one can locally use p as local coordinates, 
and make a symplectic transformation of the form (p, x) ^ {p-\- [{H"{p)) ^x, x). The tangent 
space to A is of the form {5p — [H"{p)^ ^5x, 5x) so we are reduced to replace (C, B) above 
by {C,B- {H"{p)y'c)\p^p^^), and so BC-' by SC"^ - {H")-'- By (a.6) 

(7,«) = l{sgn{BC-' - iH'r')ieo) - sgn{BC-' - (i^")" ')(-eo)) 
and as in [Ho2] this formula rewrites 

. /- \ 1^ f-H"{-eo) Id \ f-H"{eo) Id V 

This formula remains valid by continuity even if H" fails to be invertible, provided the 
matrices are non singular, i.e. so long Ty(^±^^^A is tranverse to the fiber and zero-section of 
T*M. 

B. The oscillating principal symbol and invariant densities. 

Here we review the construction of the principal symbol of asymptotic cigenf unctions 
of H, microlocalized on the invariant lagrangian manifold A, parametrized (globally) by 
(p, x) = {P{lp)^ X((/?)), Lp G T'^. For a comprehensive treatment, in the framework of Fourier 
Integral Operators calculus, we refer for instance the reader to the first Chapter of the book 
[Iv] by V.Ivrii which, together with the Lecture Notes [BaWe] by S.Bates and A.Weinstein, 
we took as a first guide to this exposition. See also [DoZh]. 

So let H{x, hDx) be the semi-classical quantization of H{j), x), i : A. ^ T*M a lagrangian 
immersion in H{p, x) — E ; we look for an approximate solution mod 0{h'^) of H{x, hDj^)u = 
E{h)u with E{h) = E + 0{h). For simplicity, we shall explicit the computations by assuming 
H{j>,x) = p^ + V{x) so that H{x,hDx) = —h^A + V{x) is a Schrodinger operator, but 
everything carries over to more general /i-PDO, as in the case of linear water waves with 
dispersion relation given by (4.10). For the moment however, we shall assume that H has no 
sub-principal symbol. 

Near xq ^ C, we try a solution of simple WKB type, of the form u{x, h) = ao{x,Ej h) 
where solves Hamilton-Jacobi equation H {x (t){x)) — E, i.e. (p^x) — pdx, or 
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d(j) = i*(pdx). This can be done so long A is projectable. To determine a{x, h) = ao{x, E) + 
ha-i {x, E) -\- ■ ■ ■ we begin with the first transport equation 

(6.2) J^0;.(x)a^,ao(x,E) + ^A0(x)ao(x,E) = 

3 

which expresses that aQ'Kx*{XH) (the x-projection of the Hamiltonian vector field of H) is 
divergence-free for the canonical density \dx\ on R^, but which we write more invariantly as 

(6.3) £x«u(ao(a;,£;)|dx|i/2)=0 

Here Cxh\k denotes the Lie derivative acting on half-densities. As Xh is tangent to A, and 
the Lie derivative invariant under diffeomorphism, this equation is satisfied iff the pull-back of 
ao(x, E)\dx\^^'^ via the projection tt is invariant under the flow of Xh- When the hamiltonian 
flow is quasi-periodic on A, and y = (p, x) = {P^ {(f) , X^ {(f)) it is readily seen that the 
solutions of (b.3) are of the form 

(6.4) ao(a;,E) = Const. Idet^^r^/^ 

d(f 

Next we analyze the situation near a focal point y, and begin with the other extreme, i.e. 
rank dTTp{y) = d. Near y, A has a simple fold and the caustic is a line ; if d = 2, then 
y ^ J-Q where Tq was defined in (a.2). We make a /i- Fourier transform, which takes formally 
H{x,hDj;) = -h'^A + V{x) to H{p,-hDp) = + V{-hDp). The action is now V(p) = 
—xdp, i.e. dip = L*{—xdp), and solves Hamilton-Jacobi equation -|- V{—Vip{p)) = E. 
So we look for u{x,h) of the form 

(6.5) I{b, iIj){x, h) = (27r/i)-'^/2 J g*(a;p+V(p))//^^(p^ e, h)dp 

and formally {H{x,hD^) - E)u{x,h) = J e^(^p+V^(p))/ft(p2 ^ V{x))b{p, E, h)dp. Using Ha- 
milton-Jacobi equation, Taylor formula gives 

(6.6) p' + V{x) -E = J2 + dpMP)) 

3 

where \j (a;, p, E) are smooth functions near x-\-dpil){p) — 0, and because we work near a simple 
fold, we can assume, possibly after a rearrangement of coordinates, that Xi\xj^q^^(^p)^q ^ 0. 
We look for a symbol h{p, E, h) = bo{p, E) + hbi{p, E) + ■ ■ ■ such that there are amplitudes 

a^^^ {x, p, E, h) = ^ (x, p, E) + htti ^ {x, p,E) + ■ ■ ■ verifying 

(6.7) e'^xp+i.{p))/hi^p2 ^ y^^^ _ ^^^^p^ E,h) = J2 hDp. {e'^^P+'^^P^^/''a^^\x,p, E, h)) 
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At zeroth and first order in h, we find respectively 

(p2 + V{x) - E)bo{p, E) = Y^i^j + dp.'^{p))a\^\x,p, E) 



ip' + Vix) - E)b,ip, E) = Y.^xj + d^,^{p))a!i\x,p,E) + 



The first equation (b.8) yields aQ \x,p, E) — bQ{p, E)Xj{x,p, E) and the second one can be 
solved iff 



r)h r)\ 

(6.9) J2>^j{-MiP),P)^iP) + ^{-dpm,P)Hp) = 

opj opj 



where we have dropped parameter E from the notations. Taking derivative of (b.6) with 
respect to Xk yields 

(6.10) d,^V{x) = d^u^A^^P^ E){xj + dp.i;{p)) + Xk{x,p, E) 

j 

and evaluating at a; + 5p(i/') = 0, 

(6.11) Xf,{-dp'^,p) = d,,V{-dp'^) 

We differentiate again (b.lO) with respect to x^, and evaluate at x + d^{p) — 0, so we get : 

(^•13) da;^Xi{-dp'ilj,p) + da:,Xk{-dpilj,p) = q^^q^^ {-dpij) 

At last, we take derivative of (b.lO) with respect to pi, and evaluate at k — £, x + dpi^) — 0, 
to obtain : 

= ^pfcAfc(-ap^,p) + ^aa^feAj(-9pV^,p) {p) 
so summing over k and using (b.l3) gives 

which we insert into (b.9) and using (b.ll) we get again : 
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This can be compared with (b.2) and again be written in an invariant form as 
{b.l5) Cx,^,{bo{p,E)\dp\'/')^0 

Since (xq) ^ 0, this transport equation can be solved with arbitrary initial condition at p = 
Pq. Again, since Xh is tangent to A, and the Lie derivative invariant under diffeomorphism, 
this equation holds iff the puUback of bo(p, E)\dp\^/'^ via the projection tt^ is invariant under 
the flow of Xh- (Here, given a lagrangian immersion i : A — > T*M, identifying locally T*M 
with R^^ we denote by Wp the projection on R^. ) 

Once we have determined bo{p,E), we can solve the first equation (b.8), and so (b.7) 
holds up to first order terms in h. Introducing a cutoff in (b.5) around po, and integrating 
(b.7) over p in the neighborhood of po, we find {H{x, hD^) — E)u{x, E, h) = 0{h'^). 

Now we pass to the general situation where rank dTZxiyo) < d — 1 and d7rp{yo) < d — 1. 
Then A can have a more complicated singularity, e.g. a Whitney crease at yo e J^o when 
d = 2, or a critical point of the moment map (which we have excluded in our special problem, 
but occurs when looking for instance at the critical energy of an harmonic oscillator, see 

e.g. [KaRo] , or more general critical points [Vu] , [CdWu] , ) From now on we drop for 

simplicity energy parameter E, and remember that everything could depend on one or several 
parameters. 

Recall from the discussion before Example a. 6 that A can be parametrized nar j/o by a 
non-degenerate phase function S(p, x) = {p, x) + H{p). More generally A can be parametrized 
by a Morse function. By this we mean a smooth function S{x,9), defined near {xo,do) £ 
R*^ xR^ (for some £ which is not necessarily the number d of phase variables, but we can choose 
i < d) with Pq = dS{xo, $q), d(^x,e)S{xQ, Oq) ^ 0, and dS'^^ • • • , dS'^^ linearly independant on 
the critical set Cs = {{x,9) e R*^ x R^ : dgS = 0}, such that the lagrangian immersions 
i : A — > T*M and ig '-Cs ^ T*M, (x, 6) i— > (x, dxS{x, 6)) have the same image near {xo,po), 
which we identify naturally with A. We set E5 = {{x,dxS{x,9)) : {x,9) e Cs} C A. Such a 
phase S{x, 6) satisfies det q-^S ^ 0, and we set cZ^^ g^S{x, 6) = S"{x, 6) (all derivatives). For 
an amplitude with compact support a{x, 6, h), we define as in (b.5) the lagrangian distribution 

(6.17) /(a, S){x, h) = (27r/i)-^/2 j ^iS{x,e)/h^^^^ ^^^^ 

and consider the /i-Fourier transformation 

J^hl{a, S){^, h) = (27r/i)-('^+^^/2 j j e'iSix,e)-xO/h^i^^^ f^^^^^g 

For ^ near po we may expand !FhI{o-i S){^,h) by stationary phase, since the critical point 
(a;(^), ^(C)), with {x{^q)., ^(Co)) = {xq, Oq) is non degenerate. The critical value equals 

(^•18) ^s{0 = s{x{0,o{0)-<0^ 
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and we find : 

(6.19) J^hI{a,S){^,h) = e''^^^^^/^\detS''\-^/^e'''^'^''^"^/%ao{x{0,0{0) + 0{h)) 

It is also well known that J^/i/(a, h) does only depends on A near {xo,po) and not on 
the particular choice of S realizing the embedding t : A ^ T*M. More precisely, if S{x, 9) 
is another phase function as above, such that po = dj:S{xQ,6Q), there exists an amplitude 
a{x,9,h) such that formally (i.e. microlocally near {xo,Po), ) I{a, S){x, h) = I{a, S){x,h). 
The convenient framework to express this invariance is described in term of the frequency set, 
see [GuSt] , [Iv] , . . . but we can (roughly) understand this by thinking of an equality between 
asymptotic sums in /i, or even at the first order in h. See e.g. [BaWe, Theorem 4.18] for a 
proof of equivalence among the Morse family ; from this follows actually Darboux-Weinstein 
Theorem, that we recalled in Sect.l. We also call S{x, 9) a non-degenerate phase function in 
the sense of Hormander. 

Choose ^ = (^1, ■ ■ ■ , ^d) as local coordinates on A and define a half density in the local 
charts (C5, is) near {xq,pq) by a/5s = | dei S"\~^/'^\d^\^^'^ . We notice that if /(a, S){x, h) = 
I{a, S){x, h) microlocally near (a;o,po), then 

(6.22) e*^s^^^/ W(e), m^l^ = e^*-(«)/^oe*'^^^^°^"-^^^^")/^ 

Let As = i^siCs) C A be the image of the chart {Cs, ls)- The (oscillating) principal symbol 
of u in A5 is then defined as 

(6.23) e**«(^)/'^ylo(e) = e'^^^^^/^e''''^''^"/^ao{xiO,diO)V^ 

To determine the principal symbol ao{x{^), 9{^)), we proceed essentially as above, computing 
H{x,hDj;) u{x,h) by inverse Fourier transform. So if the lagrangian distribution u has 
(oscillating) principal symbol e**^*^^^/'^^o(0» ^^en H{x, hDx))u{x, h), which can be computed 
as in (b.l4), has (oscillating) principal symbol e^^'^'^^^/'^ivCxH |a^o(0- (^^^ [DuHo] for the 
case of a classical PDO, and also [Iv] that takes care of specific details relative to the semi- 
classical framework. ) To have H{x^ hDx)u{x, h) = mod 0{h?), it suffices to solve 

(&.25) e'^'^^^^/^-Cx^uMOm^'^ = 

and still in the case where A is an integral manifold of H, this equation is satisfied whenever 
the puUback of the half-density ^o(0 ^^e chart (C5, is) is invariant under the fiow of Xh- 
Let us summarize our review in the : 

Theorem b.l: Let A as above be an integral manifold of H{p, x) in energy surface H{j), x) = 
E. Assume H{p,x) = p^ + V{x) so that H{x,hDx) = -h'^A + V{x), or else H{x,hDx,h) 
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is the h-Weyl quantization of a symbol of the form H{x,p,h) = H{p,x) + h'^H{p,x) (i.e. 
without subprincipal symbol Hi(p, x). ) Then near any j/o = ipoj xq) G A there is a local chart 
(C'S) i^s) and a half density as above such that the lagrangian distribution u in (b.l7) with 
(oscillating) principal part e**^(^)/''^o(0 given by (b.23) solves {H{x, hDa;, h) - E)u{x, h) = 

In particular, let y e A be such that rank d7r\A{y) = k < d — 1 is locally constant, 
then because A is lagrangian, there is a partition of variables (possibly after renumerotation 
of coordinates) x = {xi,X2) G R'^ x R''"'^, and p = (^1,^2) e R'^ x R'^~'=, such that rank 
dw = d, where tt : A — > R'', {x,p) h- > [xi,p2)- Then it is well known that u{x, h) in (b.l7) can 
be rewritten in the same form with ^ replaced by p2 and 5'o<,~-'^o7f (0^1,^2) = T{x\,p2)-\-{x2iP2)- 
In particular, near y, 

dT dT 
(6.27) A = At = {{x,p) : X2 = -^(3^1,^2), Pi = ^(0^1,^2)} 

and the number I of phase variables can be reduced to d — k. We interprete T{xi,p2) as a 
partial Legendre transformation of (p{x) or ?/'(p). See [Iv, Prop. 1.2. 5], [DoZh], . . .If we content 
computing Maslov indices, one could restrain to such phase functions as we shall explain 
below. 

Now we want to extend the previous procedure to more general /i-PDO's, in particular 
for (Weyl quantization of) H{x, hDr^, h) = Ho{x, hDj;) + hHi^x^hD^) + ■ ■ ■, in this case 
aH{x,p) = Hi{x,p) is called the sub-principal symbol of H. Of course, we need not assume 
the particular form — /i^A + V{x) for Hq. When H{x, hDx, h) is the Weyl quantization of 
symbol H{x,p) = HQ{Xjp)+hHi{x,p)+- ■ ■, then cr//(x,p) = Hi{x,p), and there are formulas 
relating subprincipal symbols for different quantizations, see e.g. [Iv]. Note that the sub- 
principal symbol, together with the principal symbol, are invariantly defined. 

Now if is a lagrangian distribution as above, then in the chart (Cg, is), the distribution 
H{x, hDxi h)u{x, h) has (oscillating) principal symbol 



^*^(^^/'^(7/:x,j. + ^h(x(0,p(0)^o(0 

L 



and again equation 



1 



{h.2Q) {-C^ I +a^^(x(0,MO))^o(OMei'/' = 



can be solved along the integral curves of |a- Where A is projectable on M, we can take 
X as local coordinates on A, and this equation takes the form 



2i ^ dpj dxj 



-Al{x) + {anlAjAlix) = 
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but near a simple fold, it can be computed as in (b.l4), and more generally using the procedure 
leading to Theorem b.l. This gives on A : 

(6.30) Ao{^) = Const.\det^\-^/^exp[-i J aHipiO<0)ds] 

Assume now that A is a torus with quasi-periodic flow of frequencies uj satisfying some 
diophantine condition we give (b.30) a simpler form. Namely expand anipi^, x{^))o[u!s + (fi) 
as a Fourier series, and integrate term by term which gives J* ch(p(C)) ^iO) ° i'^^ + (p)ds = 
{aH)t + G{(p), where G is a smooth periodic function on the torus, satisfying {G) = 0, and 
{(Jh), again, is the spatial average of an over A. Of course, this is also the time average over 
the almost periodic orbit. In (b.30) we can factor out the periodic part, which is the same 
as the solution of (b.25), multiplied by e~'^^^'^\ This will give the invariant density (since its 
Liouville class is zero. ) The other factor, e~*^^*^* is multivalued. So far we rewrite (b.30) as 

(6.31) ^o(^(^* + ^)) = Const. I det tt^I"^/^ exp[-iG{ujt + (p)] exp[-i{aH)t] 

This expression is computed along the integral curves of Xh^ ; of course it assumes a partic- 
ularly simple form when an) = 0. 

C. Maslov quantization condition for a lagragian immersion of the torus. 

At last we investigate the global properties of the lagrangian distribution u{x,h) we have 
so far constructed locally. Choose ^ = (^i, ■ ■ ■ , ^d) as local coordinates on A and define a half 
density in the local charts (Cs^t^s) near (po^xq) by ds = \det S"\~^\d$,\. We notice that if 
I {a, S){x, h) — I {a, S){x, h) microlocally near (po: ^o)-, then 

(c.l) e'^s^^^^^'ao^ = g»*.(0//.^^^g^.(sgn5"-sgn5")/4 

Let As = is{Gs) be the image of the chart {Cs,is)- We denote by fli/2 the bundle of half- 
densities on A with transition functions expz7r(sgn5'^f^ g,^ — sgnS'^^, ^-^ ^■))/4 for change 
of coordinates a; i— > a;' in A5. We denote also by L Maslov bundle on A with transition 
functions expz7r(sgn5'" — sgn5'")/4 for changes of function 5 in A5 fl A~. With I{a,S) we 
associate the section of f2i/2 (8) L of the form 

(c.2) e'^^^^^/^AoiO = e'^''^^^/^e''''^''^"/''ao{x{O,0{O)VS^ 

by a partition of unity subordinated to the covering of A by the local charts A5. 

Maslov canonical operator is defined as the inverse Fourier transform which to any section 
Qi^siO/h A(^^^fi'^ of ^1/2 (8) L as in (c.2), assigns a lagrangian distribution u{x,h). We call 
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e^^3^^y^A(^,h) the (oscillating) principal symbol of u. The interesting case is when A is 
a compact manifold. Namely, if the open sets A5 cover A, then the functions ^^^(0 = 
^siO - ^5(0 = S{xo, Oo) - S{xo, e'o), for isixo, $0) = i~{xo, O'q) = {xo,po) e A5 n A~, are 
necessarily constant on A5 n A~, defining an element $ e /f^(A;R). Actually, by De Rham 
isomorphism, $ L*{pdx). In the same way, if the covering A5 is sufficiently smooth, the 
transition functions ex ~ = e'^'^{sgnS"-agnS")/4: ^^.^ constant on A5 n A~, defining an element 
a e H^{A; Z), and a = ji{L^ L2) is Maslov class of (A, l). Of course, these transition functions 
are nothing but those we computed in Part B. Maslov quantization condition expresses the 
fact that 

- ^5(0)/^ = 7r(sgn5" - sgn,S")/4mod27r on As n A~ 
which we write as a condition on ^ L : 

(c.3) $/(27r/i) + a/4 = Omodif^(A;Z) 

The values of h satisfying (c.3) are generally of the form |- = an + 6, n G Z, and Maslov 
index a = (2, ■ ■ ■ , 2) for a torus. Global topological properties of 2-d tori are considered in 
[0-de-AlHa]. In particular, we can have a = (2, 0) in 2-d. 

At last, we point out how to correct Maslov quantization condition when H has a sub- 
principal symbol, and the frequencies of the torus are diophantine. In (b.31) e"*^'^^'** is 
multivalued, and contributes to Maslov condition by a constant term. So in case of a sub- 
principal symbol, we need to modify (c.3) as : 

(c.4) $/(27r/i) + a/4 + {aH){l, •■■,!) = Omodiy^(A; Z) 

This situation is of course well-known, and related to the sub-principal form for an integrable 
system (see [Vu] and references therein. ) 
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